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ABSTRACT 


This is the final report of a four-year numerical ocean circulation modeling prograin for the 
Gulf of Mexico. The aim of the program was to progressively upgrade an existing model of the 
Gulf of Mexico to provide the most realistic long-term simulation possible of Gulf circulation. The 
NOARL/JAYCOR multi-layer hydrodynamic and thermodynamic primitive equation circulation 
models of the Gulf of Mexico on a 0.2 degree and a 0.1 degree grid were used in the program. 
They were forced by constant inflow through the Yucatan Straits compensated by outflow through 
the Florida Straits, and by wind stresses from the Navy Corrected Geostrophic Wind data set for 
the Gulf of Mexico. The most realistic simulation produced by the program was a blending of one- 
layer and two-layer experiments, augmented by a perturbation analysis to extract surface currents 
from the vertically averaged ocean model velocities. Actual wind stresses from 1967 to 1982 were 
used to force the simulation, but the resulting currents are intended to be representative of Gulf 
circulation rather than a hindcast of the actual state of the Gulf over this time period. Wind stress, 
surface current, and geostrophic surface current fields from this simulation, sampled every three 
days for 10 years, have been delivered to MMS. Current fields at 100, 300, 750, and 1,600 meter 
depths, sampled every six days for the same 10-year period, have also been delivered. As have 
detailed vertical current profiles at the locations of 16 moored buoys used in Years | through 5 of 
MMS's recently concluded Gulf of Mexico Physical Oceanography Program. 


iv 


TABLE OF CONTENTS 


SECTION PAGE 
Disclaimer ............cccccccecccecccccccccccccceccccccececcccecccccecccccccecsees i 
PO oo n.ce'vacceevnnscnoxcneevescsececnsesese0esceuecesccevesessavswssutacresons iv 
Table of Comtent...............cccccccccccccccccccecccee, ccceeccccceceeceseess Vv 
Be OE SIN oi sosiccasecnnccnesscecssnccsansnccapscspansecscesseveuesswausiuseses vi 
Table of Model Notation...................cccccccccccccccccccecececcecceeeeeess xi 
I INTRODUCTION ............ccccccccccccccccccccccccccccccccccccccccccceceeces l 
Ul SUMMARY OF YEARS 1 AND 2 .......... 0... ccc ccccccccccccccccccccceecers 2 
Ol EMMA.  Wceenncncsianconsncdcnvccvansasnnecedquscecdnunscecsessanwasssacauacsiivenis 15 
{ YEAR  4i oii cccc cece cccccccccccccccccccccecccccccccecccceccecccccecccccccecccuceees 27 
V EP eee tO PEGI osnecaecdnndscccscsueeceseeveccedusascscssesscbdecbsunaoxecsersi 82 
REFERENCES. ................ccccccccccccccccccccccccccccccccccccccccccccccesess 212 214 
REPORT DOCUMENTATION PAGE.................cccccccccccecceeeeeees 213 


14 


15 


LIST OF FIGURES 


a. Instantaneous view of the interface deviation in a two-layer simulation 

of the Gulf of Mexico driven from rest to statistical equilibrium solely 
inflow through the Yucatan Straits (Experiment 9). b. Depth of the 
degree isothermal surface, 4-18 August 1966 (Alaminos cruise 66- 

PPR, TU COG FE Fae iceiceneceetarsccnnsansensoenecsesernissencerscsnein 


Instantaneous view of upper-layer averaged velocities from Experiment 
Ge GU GTI TR Fi i Oieccecccenenessssenensnncensncsnescccetssiosenstecssenens 


Instantaneous view of upper-layer averaged velocities from Experiment 
I TI le eitccenccetcssncesnnsccwenssnsnsicsesessmtensonssestern 


Instantaneous view of upper-layer averaged velocities from Experiment 
Oe IL ERR i iieccteeentcnnnsensnsemnenstiseesseie nnn 


Instantaneous view of upper-layer averaged velocities from Experiment 
a _, & - | | SUR IrEEeyerlinnnyrnennnr ery eT ena T TENNENT TTT 


Instantaneous view of upper-layer averaged velocities from Experiment 


201/16.0 on model Days (a) 2880 and (b) 2970.............ccccccceseeseees 
Instantaneous view of upper-layer averaged velocities from Experiment 
201/16.0 on model Days (a) 3060 and (b) 3150.................cccccceeees 
Instantaneous view of upper-layer averaged velocities from Experiment 
201/16.0 on model Days (a) 3240 and (b) 3330..0...........ccccccceeeeees 
Instantaneous view of upper-layer averaged velocities from Experiment 
201/16.0 on model Days (a) 3420 and (b) 2510...............ccccccccesseee 
Instantaneous view of upper-layer averaged velocities from Experiment 
201/16.0 on model Days (a) 3600 and (b) 3690..................cccceee 


Kinetic energy spectra at 40 m for (solid) Experiment 201/13.0 at 
Location Z ou Cashed) Mooring E sensaseendeennsndensaneteuneesenbehenntes 


Layer depths for a two-dimensional two-layer hydrodynamic model 
with full-scale bottom topography that uses Flux Corrected Transport to 
allow the layer interface to ‘intersect’ the topography....................... 


Layer on (a) Day 0 and (b) Day 3 for a two-dimensional two- 
ba ynamic simulation that uses Flux Corrected Transport to 
ow the layer interface to ‘intersect’ the topography....................... 


Layer on (a) Day 6 and (b) Day 12 for a two-dimensional two- 
a ha simulation that uses Flux Corrected Transport to 
the layer interface to ‘intersect’ the topography....................... 


Sea surface temperature and the corresponding calculated upper-layer 
density fr January from a monthly climatology based on historical ship 


vi 


PAGE 


17 


18 


19 


24 


26 


27 


31 


32 


LIST OF FIGURES (CONT'D.) 


Instantaneous (a) free surface deviation and (b) upper-layer density from 


a thermodynamic model of the Gulf with no wind forcing................. 


Instantaneous (a) free surface deviation and (b) upper-layer density from 


the same simulation as Figure 16..............ccsscsccesessscsseeseeeseceecees 


Instantaneous view of the (a) free surface deviation and (b) upper-layer 
density deviation from Experiment 212/83.1 on wind Day 093/1971 .... 


Instantaneous view of the (a) free surface deviation and (b) upper-layer 
density deviation from Experiment 212/83.1 on wind Day 183/1971 .... 


Instantaneous view of the (a) free surface deviation and (b) upper-layer 
density deviation from Experiment 212/83.1 on wind Day 273/1971 .... 


Se oe kee eatin tanacate alee 
density deviation from Experiment 212/83.1 on wind Day 1972 


Instantaneous view of the (a) free surface deviation and (b) upper-layer 
density deviation from Experiment 212/83.1 on wind Day 094/1972 .... 


Instantaneous view of the (a) free surface deviation and (b) upper-layer 
density deviation from Experiment 212/83.1 on wind Day 184/1972.... 


Instantaneous view of the (a) free surface deviation and (b) -layer 
density deviation from Experiment 212/83.1 on wind Day 274/1972.... 


Instantaneous view of the (a) free surface deviation and (b) -layer 
density deviation from Experiment 212/83.1 on wind Day 1973. ... 


Instantaneous view of the free surface deviation from Experiment 212/ 
83.1 on wind Days (a) 005/1974 and (b) 095/1974................cccce 


Instantaneous view of the free surface deviation from Experiment 212/ 
83.1 on wind Days (a) 185/1974 and (b) 275/1974.............ccccccccees 


Annual climatology from 10 years of Experiment 212/83.1 (a) interface 
deviation and (b) upper-layer density deviation......................00000 


Monthly climatology for January from 10 years of Experiment 212/83. 1 
(a) interface deviation and (b) upper-layer density deviation............... 


Monthly climatology for April from 10 years of Experiment 212/83. 1 
(a) interface deviation and (b) upper-layer density deviation............... 


Monthly climatology for July from 10 years of iment 212/83.1 (a) 
interface deviation and (b) upper-layer density ___ 


Monthly climatology for October from 10 years of Experiment 212/83.1 
(a) interface deviation and (b) upper-layer density deviation............... 


vil 


PAGE 


25 
26 
32 
33 
34 
35 
36 
37 
38 
39 
40 
41 
43 
44 
45 
46 


47 


33 


35 


37 


38 


39 


LIST OF FIGURES (CONT'D.) 


Upper-layer current vectors and bottom topography contours for the 
Texas shelf on wind Day 001/1969 for (a) one-layer Experiment 212/ 


80.0 and (b) two-layer Experiment 212/83.1 .............cccsceeeeeeeeeeeees 


Current vectors and bottom topography contours for the Texas 
on wind Day 001/1969 for (a) "geostrophic” surface currents and 


(b) matens IE crecracsn cca ncesscesnssnnnennenncnsninnentsasnansereaencees 
current vectors and bottom topography contours for the Texas 

on wind Day 091/1969 for (a) ” ic” surface currents and 

CD ID i rrsescsncenesesensnccececcesenenscctannecosescensenseessses 
current vectors and bottom topography contours for the Texas 

on wind Day 181/1969 for (a) " ic” surface currents and 

OD I GI IIE ivcceccetcnecierennsecenecceceeensnensvccsencscsosscittecsitecs 
current vectors and bottom topography contours for the Texas 

on wind Day 271/1969 for (a) ” ic” surface currents and 
i aectencctecntccecntessentaecensesnssasesttsatntvonenseccens 
Current vectors and bottom contours for the Texas 

on wind Day 002/1970 for (a) " ic” surface currents and 

CD ca neetsectetcnneicsicensnesesenceuneisescessesssesagensennees 


current vectors and bottom topography contours for the Texas 
on wind Day 092/1970 for (a) “geostrophic” surface currents and 


(a 


Current vectors and bottom topography contours for the Texas 
on wind Day 001/1969 for (a) surface currents and (b) surface 


IEE CE GD Fic I cc pen ccennscseantnnsiwetensenciessinttenennscerses 


(a) Wind Stresses and (b) local wind induced currents for the Texas 


shelf on wind Day 001/1969.........c..cccssseceseesseeecesescseeseeseeeeeees 


(a) Wind Stresses and (b) local wind induced currents from the 3.5% 


rule for the Texas shelf on wind Day 001/1969 .............c.ccceceeeeeeees 


current vectors and bottom topography contours for the Texas 
on wind Day Saat nay tO SENS CURE GRE Crp surface 


IS NE GD Bis Fe ED ccc cesccnedcecensevessnnenceceseussenseeeisonnseess 


(a) Wind Stresses and (b) local wind induced carrents for the Texas 


Ee Ok Ce By Gal re ee cewecccnsecnsescconnnsescsssnecscenssesecesess 

(a) Wind Stresses and (b) local wind induced currents from the 3.5% 

tule for the Texas shelf on wind Day 091/1969 ..............ccccceeeeeeeees 

Fide ¢ wind Day O01/1969 for tyes Experiment 
on ( i 

212/80.0 and (b) two-layer Experiment 212/83.1 .................0..cce0e0e 


vill 


49 


51 


52 


53 


61 


61 


LIST OF FIGURES (CONT'D.) 


current vectors and bottom topography contours for the Florida 
on wind Day 001/1969 for (a) “geostrophic” surface currents and 
Oy I i carcass rcccrtncnensestensneresictestermncces 


current vectors and bottom contours for the Florida 

on wind Day 091/1969 for (a) ” ic” surface currents and 

(b) Surface CUFTENIS................ccccceccceecececeececsecececesceeeeeeeeeeneees 

current vectors and bottom contours for the Florida 

on wind Day 181/1969 for (a) ” ic” surface currents and 

COD CE I isccccessesseensreerestetteeeemecnessmnnes 

current vectors and bottom contours for the Florida 

on wind Day 27 1/1969 for (a) ” ic” surface currents and 

Te I I sacoccccncccccstccncrennsevsneccnsnssuncesiecowsenseninnnees 

current vectors and bottom t contours for the Florida 

on wind Day 002/1970 for (a) ” ic” surface currents and 

(b) Surface CUPTEMIS................ccccccccccccececseccecesececscceseecececeeees 

current vectors and bottom contours for the Florida 

on wind Day 092/1970 for (a) ” ic” surface currents and 

Ce I I icemscncesaceneasccsnseserssnesenseteceeseseesseens 
Annual surface current climatology (a) from the 10-year merged simu- 

lation and (b) from historical observations ...................sccseeeeeeeeeees 


erage yy Aaenenat ri yr-toe dasdey wake dept seo x are ml 
lation and (b) “geostrophic” climatology from the merged simulation . . 


Winter surface current climatology (a) from the 10-year merged simu- 


lation and (b) from historical observations ..................cccccccseeeeeees 
Spring surface current climatology (a) from the 10-year merged simu- 
ams alae... seasewsuueeneneouuas 
Summer surface current climatology (a) from the 10- merged simu- 
tation and (0d foour tlstarical chenrvations.............-......... saskesenbsaeanous’ 
Fall surface cucrent climatology (a) from the 10- merged simulation 
ood Onoene sebenebeneenes _ sadbbnsnbasaeanavennees 

001/1969 (a) surface current from the simulation 
nd current anomaly vs the winter climatology................. 
On wind 091/1969 (a) surface current from the simulation 
and (b) current anornaly vs the spring climatology................. 
On wind Day 81/1969 (a surface cure fons the merged simulation 
and (b) surface current anomaly vs the summer. .................0.ceeeeeees 


63 


78 


73-193 


LIST OF FIGURES (CONT'D.) 


On wind Day 271/1969 (a) surface current from the merged simulation 
and (b) surface current anomaly vs the fall climatology .................... 
Map of the Gulf of Mexico showing the locations of current meter 
arrays from the Gulf of Mexico Physical Oceanography Program........ 
Velocity component spectra for 10 years from the simulation at location 
a ssckeenbensnsgbeeauecsssenenseess 
Vv component spectra for 362 days from the observations at loca- 
tion (a) u component at 60 m and (b) v component at 60 m......... 
Velocity component spectra for one year from the simulation at location 
MCD (a) u component and (b) v COmpomentt................0cceeeceeeeeeeees 
Morea fh perp son spectra for 1,090 days from the observations at 
(a) u component at 100 m and (b) v component at 100 m. 
Velocity component spectra for 10 years from the simulation at location 
MOE (a) u component and (b) v COMpomnent ...............0..ecceeeeeeeeeees 
Velocity component spectra for 10 years from the simulation at location 
MOA (a) u component and (b) v compomentl.................00c.cccceeeeeees 
Velocity component spectra for 10 years from the simulation at location 
MOC (a) u component and (b) v Compoment.................ccscceeeeeeeeees 
Velocity component spectra for 10 years from the simulation at location 
MFF (a) u component and (b) v COMpoMentt ................cccceeeeeeesreees 
Velocity component spectra for 10 years from the simulation at location 
MGG (a) u component and (b) v COMPOMERE ..............cccceceeeeeeeeeees 
Surface currents from the merged simulation, three for a full 
Loop Current eddy cycle. The wind day for cach Figure is geen by the 
"DATE" in the top right-hand corner of each plot ............0.ccccceeeeeees 


91 


Ll 


L2+L1 


L2+L1+P 


L2+L1+3.5% 


P only 
3.5% only 
FCT2D 


FCT3D 


TABLE OF MODEL NOTATION 


P.esults from a one-layer model. Currents are vertically averaged over the layer 
depth, i.e., are barotropic. The currents over the continental shelf can be taken 
as “geustrophic” surface currents. 

Results from a two-layer model. Currents in each layer are vertically averaged 
over the layer depth. The upper-layer currents can be taken as “geostrophic” 
surface currents. 

Merged L1 and L2 results using L1 in shallow water and L2 in deep water. 
Currents are still vertically averaged over layer depths. The upper-layer 
currents can be taken as "geostrophic” surface currents. 

L2+L1 with the addition of a perturbation analysis that converts the layer 
averaged currents into a full vertical profile. The actual surface currents and the 
currents at fixed depths can be obtained. 


L2+L1 with MMS's "3.5% of wind speed” rule. An alternative to L2+L1+P, 
but for surface currents only. 


(L2+L1+P) - (L2+L1) at the surface. 

(L2+L1+3.5%) - (L2+L1) at the surface. 

A two-dimensional (x-z) version of the ls yered model that used the technique of 
“Flux Corrected Transports” to allow layer interfaces to effectively intersect the 
bottom topography. 

A three-dimensional version of FCT2D. This model was never completed. 


Clima. logical currents from historical ship drift data. 


I. INTRODUCTION 


This report gives details of the third and fourth years of a tour-year numerical ocean 
circulation modeling program for the Gulf of Mexico. The aim of the program was to 
progressively upgrade, in modest increments, an existing numerical ocean circulation model of the 
Gulf so that the final model has a horizontal resolution of about 10 km and a vertical resolution 
approaching 1-10 m in the mixed layer, 10 m at the thermocline, and 100 m in the deep water. 
Throughout the four-year pesiod, the validity of the upgraded model was continuously tested and 
velocity field time series were delivered periodically based on the most realistic simuiation of Gulf 

Experiments in the first year were with the existing Naval Oceanographic and Atmospheric 
Research Laboratory* and JAYCOR (NOARL/JAYCOR) two-layer hydrodynamic primitive 
equation ocean circulation model of the Gulf on a 0.2 degree grid. They concentrated on correctly 
specifying the coastline and bottom topography for maximum realism in circulation simulation and 
on how best to include wind forcing. 

Experiments in the second year were also with iiie existing NOARL/JAYCOR 
hydrodynamic model, but most were on a 0.1 degree, rather than 0.2 degree grid. ‘Wind-forced- 
only, port-forced-only, and wind-plus-port-forced simulations were generated. The major 
limitation of layered models is that bathymetry is confined to the lowest layer, so multilayer models 
are not generally realistic in shallow water. Detailed comparisons between the model results and 
observations on the Florida Shelf demonstrated that this limitation was severe enough to warrant a 
change of approach. 

An attempt was made in the third year to modify tne layered model to allow bathymetry to 
be in any layer, and hence, more realistically represent shelf regions. A two-dimensional prototype 
was successfully demonstrated, but the three-dimensional version was not functional by the end of 
the third year. 

Rather than continue the development of this new model design, the Year 4 effort 
concentrated on producing the best possible surface currents with the existing thermodynamic 
model. This involved running a one-layer version of the model to obtain realistic results over the 
shelf and blending them with two-layer model results that were more realistic in deep water. A 


* NOARL was previously known as the Naval Ocean Research and Development Activity (NORDA). 


perturbation analysis technique was then used to obtain vertical profiles from the blended layer 
averaged fields. 

Simulated surface currents sampled every three days for more than 10 years were delivered 
to MMS from the blended and perturbed fields, which represents the best simulation available from 
the entire four-year effort. Corresponding currents at selected fixed depths sampled every six days 
were also delivered, as were detailed vertical profiles of currents sampled every three days at 16 
locations where MMS has collected current meter time series in the Gulf. 


IL. SUMMARY OF YEARS I AND 2 


The hydrodynamic ocean model, used for Years 1 and 2, contained many innovative 
features and the reader is referred to its detailed discussion in Hurlburt and Thompson (1980). In 
particular, Section 2 (pp.1613-1614) gives the model equations and Appendix B (pp. 1647-1650) 
describes the numerical design of the model. Since 1980, the code has been completely rewritten 
by JAYCOR and © «= <p. vility to handle general basin geometry had been added, but this does not 
affect the descri. on in any major way. Wind forcing is treated identically to interfacial and 
bottom stress terms, ..c., wind stress appears directly as an additive term in the momentum 
equation [p.1614 of Huriburt and Thompson (1980)]. 


In terms of ‘realism’ Experiment 9 (Wallcraft, 1984) was the most successful Gulf of 
Mexico numerical simulation prior to the start of this project. The model was driven from rest to 
statistical equilibrium solely by a steady inflow through the Yucatan Straits that was compensated 
by outflow through the Florida Straits. Figure 1 compares ‘instantaneous’ upper-ocean flow 
patterns just before an eddy is shed by the Loop current (a) from the numerical model and (b) from 
observations by Leipper (1970). The ability of the model to simulate observed features is clearly 
demonstrated by this comparison, which is remarkable given the simplicity of the model forcing. 
However some discrepancies remain; for example, the eddy has not penetrated as far into the Gulf 
and is more intense than that shown in the observations. Waves can be seen moving around the 
wall of the Loop in the observations, but in the model they are at the limit of resolution and 
therefore unrealistically large. Moreover in the Gulf the waves are more pronounced on the eastern 
wall of the Loop and can form strong cold intrusions that may contribute to the eddy-sheddir- 
process (Vukovitch and Maul, 1984). 


Experiment 68 (Wallcraft, 1984) is the most realistic simulation produced in Year 1. It ia. 
two layers on a 0.2 degree grid and is forced by transport through the Yucatan Strait and by winds 


Figure 1. (a) Instantaneous view of the interface deviation in a two-layer simulation of the Gulf ci 


(a) 


(>) 


Mexico driven from rest to statistical equilibrium solely by inflow through the Yucatan 
Straits (Experiment 9). The contour interval is 25 m with solid contours representing 
downward deviations. (b) Depth of the 22-degree isothermal surface, 4-18 August 1966 
(Alaminos cruise 66-A-11), from Leipper (1970). The contour interval is 25 m. 
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from the Navy Corrected Geostrophic Wind data set for the Gulf of Mexico, which is based on the 
Navy's 12-hour surface pressure analysis from 1967-1982 (Rhodes et al., 1986). Figures 2-5 
show upper-layer currents (i.e., vertically averaged currents above the thermocline) every 60-120 
days for 300 days. Vectors are only drawn at every second point (i.e., every 0.4 degrees) to 
improve readability. Figure 2 shows the furthest northward penetration of the Loop Current ever 
attained by the ocean model, which is often seen in the Gulf. After the eddy breaks off, the Loop 
Current intrudes onto the Florida Shelf and some of the flow splits off to the north for a brief time 
(Figure 3). Similar intrusions have been observed in the Gulf, but the model's inadequate 
representation of shelf topography makes it likely that the simulated currents in shallow areas (less 
than 100 m) are too energetic. A persistent anticyclonic gyre in the northwest Gulf has been a 
feature of almost all Gulf simulations performed to date. The addition of wind forcing in 
Experiment 68 has increased its average size and its effect on incoming Loop Current eddies 
(Figure 4). The presence of a gyre in this position is explainable by the northward migration of 
anticyclonic eddies along the coast of Mexico until the continental slope bends eastward and they 
can go no further, which magnifies the effect because the winds also tend to produce an 
anticyclonic gyre at the same spot. However, in the Gulf the gyre probably dissipates relatively 
rapidly against the shallow shelf area. 


Simulations in Year 2 (Wallcraft, 1986) were on a 0.1, rather than 0.2 degree grid. 
Experiment 201/16.0 is the most realistic of these and it is forced by winds and Yucatan transports 
very similar to those used in Experiment 68. Figures 6-10 show upper-layer currents every 90 
days for two Loop Current eddy shedding cycles with winds from 1968 and 1969. When 
compared with similar simulations without wind forcing, this simulation has higher maximum 
currents and significantly more variability; for example, in Loop Current eddy shedding cycle 
times. However, the anticyclonic gyre off the Texas coast is generally larger, as is expected from 
the wind stress curl patterns. 


The major limitation of layered models is that bathymetry is confined to the lowest layer, so 
multilayer models, in general, are not realistic in shallow water. In order to quantify how well the 
ocean model was performing in shelf areas, three years of hydrographic data obtained at fixed 
moorings on the Florida Shelf (as part of a 1986 MMS-sponsored Gulf of Mexico Physical 
Oceanography Study done by SAIC) was compared to 1,080 days of simulated current meter 
data at approximately the same 1/13.0. In all cases the currents were significantly stronger in 
the simulation than in the observations. Figure 11 shows a comparison of observed and simulated 


Figure 2. a view of upperdayer averaged velocities from Experiment 68 on model 
Day 3918. Vectors are only plotted at every second model grid point, i.e., every 0.4 


degrees. 
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Figure 3. Instantaneous view of upperlaycz averaged velocities from Experiment 68 on model 
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Day 4038. Vectors are only plotted at every second model grid point, i.e., every 0.4 
degrees. 


GE@STR. CURRENTS G. OF MEXICZ 0. 6&8 


MZOEL ORY = 4038 WINO ORY = 1968/054 


~e4& & AD a 


i 


! i ‘ ‘ ‘ ‘ ' i i i 4 i i i ' i i 
: : : J : | Jan > > rs : 
' , £4eravyr»v 
: rZeneR a4 <¢€¢ £7 PD? 
—— Viea & o2<Preecastt®Fagrr - 
: ee a Ye e 24M ~ 
4% 44¢a4¢ —n—eo meee ect t & & &Y x ew <lv, . 


_ 
‘ 4 2 A-One, & > F 
\Nweruag es? 
4s - 
al a at Ean 
ee ala as = ~ 
' 
' 
! L L ! ! ihe 


Core SC eESTOOSOMRGD TC KR? OMFIF ED CHV 
Seat £7 F Meee tH Geese tt Aw Set yg Vitis 


vo 44 © tee tara rivers 
<4 4f7Z ie ae ee 


aay, “at A 72 th asee~< ‘ ' 
vig <a" F 7 o » > ~ ~ e° 
eer LA - v6“ vf : 


3 et OG 


? 
wot teasing 
a> 4 
eXN 


7 
NTR LZ, ere eet ec rian oo 23 


at Laat en oem parce : . 
* jem) ah titearovsveaaat 4> 4, ‘ 
at + agi rat aaemuaueé * _ 

a 


aera ereTeriTe ST eo a 0 


MAX PL@TED VECT@R = 1.01 (M/SEC) 


Figure 4. Instantaneous view of upper-layer averaged velocities from Experiment 68 on model 
Day 4158. Vectors are only plotted at every second model grid point, i.e., every 0.4 


degrees. 
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Figure 5. Instantaneous view of upper-layer averaged velocities from Experiment 68 on model 
Day 4218. Vectors are only plotted at every second model grid point, i.e., every 0.4 


degrees. 


GEZSTR. CURRENTS G. OF MEXIC2 CO. 68 
MODEL ORY = 4218 WIND ORY = 1968/2324 


SIN ' ‘ ' La is ‘ ot i ig i t ‘ i i 


] ' ' , ' | ' see ' 
‘ id 
' <r? 77% awe 


! 

i 

' 

u pA Sf FP eOraqgeeer«< 


yee Ph Ks aces 

at XN et ewadt 
74 pn a 4 a i 
aeer"taga* 


ax*<V¥eygrr 


RHEE N 

(DEG) Wescaeen ] 
anette ae eee 

Sa 


Perea n art araav 
ee ale ap eX\&Ceeeraa 
Perens © 


L8N | L L 
98H (DEG) SOW 


MAX PLETED VECTOR = 1.84 (M/SEC) 


Figure 6. Instantaneous view of the upper-layer averaged velocities from Experiment 201/16.0 on 
model days (a) 2880 and (b) 2970. The assigned dates (344/1967 and 069/1968) 
indicate the applied wind forcing, but the experiment was not a hindcast and the ocean 
currents in the Gulf on that date might have been quite different from those shown. 
Vectors are only plotted at every second model grid point, i.e., every 0.2 degrees, and 
all velocities greater than 50 cm/sec are plotted as 50 cm/sec. 
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Figure 7. Instantaneous view of the 


yer averaged velocities from Experimen 


t 201/16.0 on 


upper-la 
model days (a) 3060 and (b) 3150. The assigned dates (159/1968 and 249/1968) 


indicate the applied wind forcing, but the experimen 


t was not a hindcast and the ocean 


currents in the Gulf on that date might have been quite different from those shown. 
Vectors are only plotted at every second model grid point, i.e., every 0.2 degrees and 
all velocities greater than 50 cm/sec are plotted as 50 cm/sec. 
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Figure 8. Instantaneous view of the upper layer averaged velocities from Experiment 201/16.0 on 
model Days (a) 3240 and (b) 3330. The assigned dates (339/1968 and 063/1969) 
indicate the applied wind forcing, but the experiment was not a hindcast and the ocean 
currents in the Gulf on that date might have been quite different from those shown. 
Vectors are only plotted at every second model grid point, i.e., every 0.2 degrees and 
all velocities greater than 50 cm/sec are plotted as 50 cm/sec. 
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Figure 9. Instantaneous view of the upper layer averaged velocities fror.. Experiment 201/16.0 on 
model Days (a) 3420 and (b) 2510. The assigned dates (153/1969 and 243/1969) 
indicate the applied wind forcing, but the experiment was not a hindcast and the ocean 
currents in the Gulf on that date might have been quite different from those shown. 
Vectors are only plotted at every second model grid point, i.e., every 0.2 degrees and 
all velocities greater than 50 cm/sec are plotted as 50 cm/sec. 
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Figure 10. Instantaneous view of the 


upper layer averaged veliccities from 


t 201/16.0 


on model Days (a) 3600 and (b) 3690. The nasigned dates (333/909 and 058/1970) 
indicate the applied wind forcing, but the experiment was not a hindcast and the ocean 
currents in the Gulf on that date might have been quite different from those shown. 
Vectors are only plotted at every second model grid point, i.e., every 0.2 degrees and 
all velocities greater than 50 cm/sec are plotted as 50 cm/sec. 
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Figure 11. Kinetic energy spectra at 40 m depth for (solid) Experiment 201/13.0 at location Z and 
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kinetic energy spectra at mooring MOF, which was located on the inner shelf in SO m of water. 
The simulation is far too energetic and still shows a strong influence from the Loop Current that is 
completely absent from the observations. 


ITI. YEAR 3 


A major thrust of the third-year program was an attempt to produce improved simulations 
in shelf areas by developing a version of the layered-ocean model that would allow layer interfaces 
to, in effect, intersect the bottom topography. When intersection occurs in a conventional layer 
model the layer thickness becomes negative, which is clearly unphysical and leads to unrealistic 
results and, if the situation persists, catastrophic failure of the model run due to undamped 
instabilities. The obvious solution of setting a minimum layer thickness at o- about zero does not 
work because: (a) it leads to loss of mass and (b) clamping the layer thickness induces dispersive 
ripples in the interface at the intersection point (i.c., we have an unresolved boundary layer). In a 
study of the effects of hurricanes in the Gulf, Cooper and Thompson (1989) used a multilayer 
model in which each layer interface intersects with a different artificial vertical step in the 
topography that acts as a solid wall boundary to that interface. This technique proved effective for 
short-term simulations of hurricane forcing with the simplifying assumption that the Gulf was at 
rest before the hurricane occurred. The technique is not generally applicable to simulations that 
include the Loop Current. In this case, the variations in each layer interface depth across the Gulf, 
which makes it impossible to design a stepped topography that guarantees that each interface will 
be everywhere and intersect the topography on the vertical face of a step at all times. The most 
promising approach to the layer intersection problem is to insure positive layer thicknesses via 
‘Flux Corrected Transport,’ a technique that was originally developed for fluid problems with 
shocks (Book et al., 1981). In this method the continuity equation is solved for layer thickness 
using two different sets of wansports; one obtained via a low-order (highly dispersive) scheme, 
guaranteed to give monotonic results; and, the other via a standard (ripple prone) high order 
scheme. The low-order scheme, used alone, would prevent layer intersection, but it very rapidly 
damps out circulation features and would not produce realistic simulations. Instead, the final layer 
thickness at each point is a linear combination of the two solutions chosen to be as close as 
possible to the high-order solution. Away from areas of layer intersection, the high-order scheme 
is used alone and the solution is identical to that without FCT, but near intersections just sufficient 
contribution from the low-order scheme is used to ensure a positive layer thickness. In other 
words, bottom topography is still confined to the lowest layer, but that layer can get very thin so 
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there is effectively no contribution from the deep layer over the shelf and no limit on how shallow 
the bottom topography can be. This method has been used with some success in a similar layered 
ocean model, both for interfaces that intersect the surface and for layers that intersect the 
topography (Bleck et al., 1983). The major problem with the method is that FCT is an inherently 
explicit scheme in contrast to the existing ocean model that treats gravity waves impiicitly (to allow 
much larger timesteps). A fully explicit model's timestep would be controlled by the external 
gravity wave speed (about 150 m/s), but the depth averaged flow can still be treated implicitly in a 
model with FCT, so the timestep depends on the internal gravity wave speed (about three m/s). 
The ocean model without FCT, on the other hand, treats both external and internal gravity waves 
implicitly and can use a timestep three-to-five times longer than "he FCT code. 


Two-dimensional (x-z) versions of a two-layer hydrodynamic model that uses FCT to 
allow layers to intersect the bottom were developed and tested on sections across the Gulf of 
Mexico (FCT2D) (Figure 12). The position of each model grid point is indicated by a vertical line 
below the topography contour. Data is only available at grid points and straight lines are used to 
connect data values in the plot. The upper-layer rest depth is 300 m over deep water, but is less 
near 98W and 82W, where the continental shelf is shallower than 300 m deep. The lower layer is 
set to be at least 10 cm thick across the entire region, so there is a lower layer over the continental 
shelf although it is too thin to be seen in the plot. Figure 12 is for two days into an experiment to 
test the ocean model with no applied forcing. The layers are in exactly the same position as at the 
initial time and the velocities are zero everywhere. This demonstrates that the model does not 

Figures 13 and 14 show only the upper 450 m of the water column for a gravity wave 
sloshing experiment, where there is no applied forcing, but the layer interface is initialized with a 
single-period cosine profile across the region. Figure 13(a) shows the initial state with about 100 
m variation in the depth of the interface from east to west (note that the lower layer is again 10 cm 
thick where the topography is shallower than the expected interface depth). Figure 13(b) is for 
Day 3 of the simulation, where the layer interface is now almost level. Figure 14(a) is for Day 6. 
The layer interface has moved up or down about 100 m at each end to reverse the profile; however, 
it is no longer exactly sinusoidal because gravity waves travel slower in shallow water than they do 
in deep water. The interface is level again between Days 9 and 10. Figure 14(b) is for Day 12, 
where the interface is again shallower to the west as it was on Day O, but the wave is almost square 


Figure 12. Layer depths for a two-dimensional two-layer hydrodynamic model with full-scale 
bottom topography that uses Flux Corrected Transport to allow the layer interface to 
‘intersect’ the topography. The figure is for Day 2 of an experiment testing the 
stability of the rest configuration in the absence of external forcing. There has been no 
change over the two days. The lower layer is 10 cm thick at all points where the 


topography appears to intrude into the upper layer. 
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Figure 13. Layer depths on Days (a) 0 and (b) 3 for a two-dimensional two-layer hydrodynamic 
simulation that uses Flux Corrected Transport to allow the layer interface to ‘intersect’ 
the topography. Only the upper 450 m of the water column is shown. 
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Figure 14. Layer depths on Days (a) 6 and (b) 12 for a two-dimensional two-layer hydrodynamic 
simulation that uses Flux Corrected Transport to allow the layer interface to ‘intersect’ 
the topography. Only the upper 450 m of the water column is shown. 


LAYER DEPTHS COM. FCT 220322 783 


Y= 26.00N Day = 6000 


3300 1 t | i tt ¢@ ¢ 4 i i i oF i 
(i 
Aiiit 
| { | 
| | 
| 
Cs | 
| 
2800 L I 1 i l t 1 t I a a | 
oew (DEG) 62W 


LAYER DEPTHS C.0.M. FCT 22082:2: 75.5 
Y= 26.00N DAY «= 12.000 


3300 a on ee ee ee ee ee ee ee ee ee ee 
A 
| if 
/ 
| 
anee Lippy py yp py yp gy yg yf 
oav (DEG) 62w 
19 


BEST COPY AVAILABLE 


and the model blows up at Day 15 as the wave ‘breaks.’ The conversion of the original wave into 
a breaking wave is to be expected, given than gravity waves travel slower in shallow water. 


The extension of FCT2D to three dimensions (FCT3D) was still under development at the 
end of Year 3. The only difficulty encountered was the practical one of the time and expense 
involved in debugging the code. There is no reason to suppose that a model of this kind could not 
work in the Gulf of Mexico; however, there is no guarantee that it could produce acceptable 
simulations. One obvious area of concern with such a model would be the accuracy of the 
simulation in the zone where the model goes from one layer to two. However, any single model 
that tried to accurately simulate on- and off-shelf flow at the same time would have problems in this 
shelf slope region. Models like FCT3D are at least one, and perhaps two, orders of magnitude less 
expensive than all other model designs that might be able to provide similar simulations. 


The second major element of the third-year program was the addition of bulk-layer 
thermodynamics to the model. The existing model was ‘hydrodynamic’ (i.e., the density in each 
layer was constant and fixed for all time) with the fixed densities and the density difference 
between layers set to realistic mean values for the Gulf of Mexico. With the addition of 
thermodynamics, the density is allowed to vary in space and evolve with time under the control of 
an equation of state and a temperature equation forced by heat fluxes (or by relaxation back to a 
density climatology) and the density of inflow water. Mixing of denser fluid from the layer below 
into each layer is also allowed and this can prevent the layer interface surfacing due to upwelling, 
i.e., the layer gets denser rather than getting thinner. In hydrodynamic models there is no 
exchange of fluid between layers. 


The vertically integrated equations of motion used in the n-layer finite depth thermodynamic 
model are, for k=1...n: 


av, 


a t(V Uti Vo+kx (% = 
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° l=1 


i=l 
+(Aa1 — 7.) / Po + AgV* ty 
Ope . . 
ny +t: Ven = 


max(0, un)“7—"* + max(0, 4-1) S 


20 


= H . 
— Bolle «(on - be) - bua + Kuen 


Pury Picphy 
Gh P 
= +V-VYy,= 
hy ( Pelle ‘ 19 ) 
— wey + — | - 6 
Wk — Wer + re ce 7,( Pe — Pr) + ens 
where 
Vi = heir 
D = actual water depth 
D = model water depth, D > D 
H, = k-th layer thickness at rest 
n-1 
H, = D(z,y)-> &, 
isi 
h, = min(h,,D) 
Pe = k-th layer (constant) density at rest 
f,» = k-th layer density climatology 
Q = surface heat flux 


The above equations are only valid provided the layer thicknesses are always positive, and 
the density fields are always stably stratified. Stable stratification is not generally a problem in the 
Gulf, and mixing can be used to prevent a layer becoming thin due to upwelling, but the model 
depth, D, must be made deeper than the actual Gulf depth, D, over the continental shelf in order to 
insure that the lowest layer interface never intersects the bottom. 


The model is always initialized from rest, with zero velocities, level layer interfaces, and 
constant density fields, i.e., at time t=0 : 


Vy. = 0 
Pr = Ph 
h, = Hp 


In the absence of external density forcing and mixing, the density in each layer remains 
constant for all time and the above equations are equivalent to those for the hydrodynamic model 
(Hurlburt and Thompson, 1980). External density forcing can take the form of atmospheric heat 
flux, Q, or relaxation back to a density climatology, f,, or the inflow through a port can be at a 
different density. Mixing is the process of entraining or detaining mass from one layer to the next. 
Mixing cannot be applied to a hydrodynamic model because it causes a change in the layer 


21 


averaged density. In a layered model, the “physics” of mixing is controlled by the calculation of 
the vertical mixing velocity at the layer interfaces, Oy. Since in this case the layer interface 
represents the thermocline depth, rather than the mixed layer depth, very simple mixing physics is 
used to effectively constrain the layer thickness to tend to stay in the range hy to hg. This form of 
mixing is primarily designed to prevent the layer interface surfacing due to upwelling, i.e., the 
layer gets denser rather than getting thinner. Mixing is also balanced, by Oy, so that there is no net 
exchange of mass between the layers over the region as a whole. 


The surface heat flux fields, Q, used to drive an ocean model must be of high quality, 
otherwise they will cause climate drift in the model by applying too much heating or cooling. 
Since there are no data sets of sufficient quality for heat fluxes in the Gulf of Mexico, the model 
was instead required to relax back to a monthly density climatology, A, based on sea surface 
temperature maps from historical ship observations. One way to look at relaxation to climatology 
is as a heat flux that always drives the model's instantaneous density field back towards a reference 
density field, i.e., 


Pifty . 19 
—— 0 — => 


where 


Q = EiSpie ae (p3 - in) 


The e-folding time for the relaxation, ph,/p;Hi Gp, was chosen to be about 360 days in 
the center of a Loop Current eddy, but only 90 days over the continental shelf. The use of the 
smaller of the layer thicknesses and the actual water depth, hy, ensures rapid relaxation over the 
continental shelf where rapid winter cooling is important. The model can be forced by sufficiently 
accurate heat flux fields; however, because of the problem of climate drift, it would probably 
always be wise to include a relaxation term that forced the model's density field to stay close to the 
climatological mean field. This relaxation term is, in effect, a parameterization of the mean heat 
flux forcing. The additional heat fluxes, applied to a model that includes such a relaxation term, 
should represent the perturbations in heat fluxes about the long-term mean. The present model is 
forced by relaxation to a monthly density climatology. The alternative, should heat fluxes of 
sufficieit quality be available, would be to relax to a mean density climatology and in addition 
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force the model with heat flux anomalies about the mean heat flux field. The mean density 
climatology would simply be the average of the existing monthly density fields. 


Figure 15 shows the sea surface and corresponding density fields for January. The density 
ranges from 25.8 sigma-t over the shelf to 24.0 sigma-t in the Yucatan Straits. The model sets the 
density of the Loop Current inflow through the straits to agree with the density climatology there 
month by month. Figure 16 shows contour maps of sea surface height and of density for an 
experiment without wind forcing and with a constant inflow transport of 20 Sv in the upper layer. 
It is for a day in December and so the shelf areas are cool (note the strong density front between the 
large anticyclonic and the smaller cyclonic eddies in the northwest Gulf). In this case, a new eddy 
is in the process of developing on the Loop Current as it extends into the Gulf. Figure 17 is for 
180 days later in the summer. The anticyclonic eddy has detached from the Loop Current and is in 
the central Gulf, which is already cooler than the Loop Current. The shelf areas are significantly 
less dense than they were during the winter. 


A detailed study of the appropriate model parameters to use for this layered model in the 
Gulf of Mexico was performed by Hurlburt and Thompson (1980 and 1982). Their results were 
for a simplified Gulf of Mexico geometry, but they carry over to the more realistic geometry used 
here. Many of the parameters are constrained to be physically realistic for the Gulf of Mexico. 
The Coriolis parameter and the beta used in the beta-plane approximation are set by the Gulf's 
location and the acceleration due to gravity is set at 9.8 meters per second for obvious reasons. 
The mean upper-layer thickness, H1, is set to 200 m in all two-layer simulations because this is the 
approximate average depth of the thermocline in the Gulf. Similarly the value of reduced gravity, 
g', is chosen based on the average contrast between the vertically averaged density above and 
below the thermocline. The allowable range is from 0.02 to 0.03 meters per second and 0.025 is 
typically used here. The reference thickness of the lower layer, H2, in two-layer simulations is 
chosen to be slightly larger than the actual rest thickness everywhere when allowing for bottom 
topography. Since the bottom topography used is never deeper than 3,650 meters, a value of 
3,450 meters is used for H2 (allowing for an H1 of 200 meters). The actual rest thickness of the 
lowest layer depends on the bottom topography. Following Hurlburt and Thompson (1980), the 
total transport though the Yucatan Straits is held fixed in all multilayer experiments at 30 Sv; 
however, the distribution of transport between layers is difficult to estimate from observations and 
different simulations have allocated as much as 26 Sv and as little as 20 Sv of the total to the upper 
layer. Allocating more transport to the upper layer increases the speed of the Loop Current; how- 
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Figure 15. Sea surface temperature and the corresponding calculated upper-layer density for 
January from a monthly climatology based on historical ship observation data. 
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Figure 16. Instantaneous (a) free surface deviation and (b) upper-layer density from a 
ic model of the Gulf with no wind forcing. Snapshots are for December 
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Figure 17. Instantaneous (a) free surface deviation and (b) upper layer density from the same 
simulation as Figure 16. Snapshots are for 180 days later than Figure 16, i.e., for a 
day in the summer. 
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ever, allocating more transport to the lower layer tends to prevent the Loop Current from intruding 
onto the Florida Shelf. The horizontal grid spacing is chosen to resolve the features of interest. A 
good rule of thumb is that 10 grid points are required across and eddy and five grid points across a 
jet in order for their dynamics to be resolved by the model. Twenty kilometers is sufficient to 
resolve the Loop Current and Loop Current eddies, but 10 kilometers is required to resolve some 
associated features that have smaller horizontal scales. The actual grid spacing is specified in 
degrees at the center of the Gulf and is either 0.2 degrees (20 km by 22 km) or 0.1 degrees (10 km 
by 11 km). The model parameterizes subgrid scale processes using Laplacian eddy viscosity and 
density diffusivity. The appropriate coefficients are primarily controlled by the horizontal grid 
spacing and less viscosity is required on a finer grid. Density diffusivity is only required in 
thermodynamic models, since hydrodynamic models have no horizontal variation in density. 
Density diffusivity has typically been set five times larger than eddy viscosity. A quadratic 
formulation has been used for bottom stress with a typical coefficient of 0.002, although some 
simulations have used 0.003. 


IV. YEAR 4 


The thermodynamic simulations started in Year 3 were continued into Year 4; however, no 
further attempt was made to develop a three-dimensional model with flux-corrected transport. 
Instead, the best possible simulation was generated using the existing multilayer thermodynamic 
model. Similar barotropic models have been used by others (O'Brien, 1971). The 
NOARL/JAYCOR-layered ocean model is designed to run with one or more layers. The one-layer 
Gulf of Mexico simulations use the same model code on the same grid over the same region as the 
two-layer simulations. The only difference is that there is just one layer and the topography does 
not need to be scaled back as it does in multilayer simulations. 


The model equations apply equally well to a one-layer barotropic model as to a two- or 
three-layer model. In fact, the original Gulf of Mexico simulations with a very early version of the 
model included barotropic experiments (Hurlburt and Thompson, 1980). The advantage of one- 
layer simulations is that the standard requirement, rather than bottom topography confined to the 
lowc st layer, is not a limitation (since the only layer is the lowest layer). The disadvantage of such 
a model is that one vertically-averaged layer is a very poor approximation to the ocean's structure 
in deep water since it cannot include a thermocline. So a one-layer model does well over the shelf, 
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but less well over deep water; and a two- or three-layer model does well over deep water, but not 


as well over the shelf. 


To obtain the “best” simulation, the results from two independent model runs (one layer 
(L1) and two layers (L2) using the same forcing) were merged (L2+L1). The barotropic fields 
were used in water shallower than 100 meters, the two-layer fields were used in water deeper than 
200 meters, and a linear combination of the two sets of fields weighted by water depth were used 
in the fairly narrow zone at the shelf break between these regions. In particular, the merged layer 
thicknesses and transports are given by: 
girth) - avi”) +(1- a)o”? 


gllitla) _ 0 for D< pio) 
: 7 ai”) for D > ni) 
where 
l for D < 100 
ars 0 for D > 200 
(200 — D)/100 otherwise 


and the merged layer thicknesses are given by: 


pliarrea) - min(h{) D) 


piirta) _ 0 for D < Al) 
A?) for D > Ai? 


The choice of 100 m as the cutoff depth, at which L1 starts merging with L2, was chosen 
based on current meter data on the Florida shelf (SAIC, 1986) that indicated that there was no 
discernable effect of the Loop Current in water shallower than this depth. The choice of 200 m as 
the cutoff depth, at which L2 starts merging with L1, was based on the average depth of L2 layer 
interface across the entire Gulf at this depth. Any value between 200 m and 500 m would have had 
very little effect on the final result since the horizontal distance between these depths is very small 
in the Gulf. 

This approach does not guarantee that currents in the blending zone will be physically 
realistic. For example, an attempt to merge Experiment 210/16.0 from Year 2 (Figures 6-10) with a 
barotropic experiment would give poor results since modeled deep water featun:s often intrude to a 
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significant degree onto the inner shelf. These features, such as an anticyclonic eddy off the Texas 
coast and the Loop Current on the Florida shelf just after an eddy is shed (Figure 7) would be 
sheared off in the merged fields. Merging the fields only makes sense if the multilayer simulation 
confines dee; water circulation features primarily to deep water. Three-layer models have 
difficulty keeping circulation off the shelf because there is insufficient connection between the 
topography confined to the lowest layer and the upper-layer circulation. Two-layer models 
typically do better than three-layer versions at confining flow to deep water, but conventional two- 
layer simulations are not good enough to allow merging of the one- and two-layer fields. A new 
two-layer version of the model was used with a modified mixing term that holds the layer interface 
close to the rest depth of 200 meters over the shelf for the entire ‘mulation. Major eddies and jets 
cause significant deformation of the layer interface, so these fea .-e cannot get over the shelf if the 
interface is held rigid. The advantages of this approach are that the minimum topography depth can 
be raised from 500-300 meters and the simulation is much improved in the shelf and shelf break 
regions. The disadvantages are that only two layers can be used (to maximize shelf realism) and 
the strong mixing over the shelf must be balanced region wide to maintain a constant-layer 
thickness that adds friction to the deep water circulation. In the Gulf, there are occasional 
intrusions of eddies onto the outer shelf and more frequent entrainment of fluid off the shelf into 
deep water. The former are excluded from the merged simulation, but the latter can occur since the 
two-layer component includes the full shelf region. A major advantage of blending two 
independent simulations is that either one can be replaced by an improved simulation at a later date. 
In this case, the shelf circulation is more likely to be upgraded in this way. In fact, each major 
shelf region (i.e., the Texas and Florida shelves) could be upgraded with different models at 
different times. 

Layered models provide current fields representing the vertical average over their 
instantaneous layer depths. Upper-layer currents are very similar to geostrophic surface currents. 
In order to obtain actual surface currents and currents at fixed depths throughout the water column, 
an analytic perturbation analysis is used (Thompson, 1974). By subtracting the vertically-averaged 
equations from the Navier-Stokes equations, diffirentiating the result with respect to the vertical, 
employing the hydrostatic relation, and assuming small Rosby number and quasi-equilibrium 
conditions, the following lowest order perturbation equation results: 


a = Oz,y,t)+a'(z,y, z, t) 


ee 


. dv 
kx f— = ~ Vp + Avs 
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wy = (U, + u,) +i(d, + v;,) 0<z7<h =—-hA, 


oe _ om 
6z3 Oz. 


and boundary conditions 
Oz | 2 Pedy 
Owe Owe; 
Bele, = OF key, "2 BRR 
Wl or, = Weil ate = J, kel...a-1l 
= Av dw. =J 
alee OO le 


The boundary conditions match both stress and velocity across each internal interface and 
linear friction is applied at the bottom, i.e., Ay w/z = Cow. 


If all the Qys and Js were known, the solution would be: 


ao 23S sa 
s s s? 
Ae = 0-8 +G./s" 
BR = a) Ae 
% 
+0, ( eet) + a - H 


G 
CO = Aue) + Fel) - Sh +m 


The boundary values Q;...Q, and J)...J, can be found by solving (numerically) an appropriate 
set of simultaneous linear equations formed by expressing A, B, and C in terms of Q and J, and 
expanding the boundary condition equations. 

The resulting currents (L2+L1+P) have the same vertical average as the layered fields 
(L2+L1), but include an Ekman component and a contribution from horizontal variations in 
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mass. The major disadvantage of this approach is the assumption that the model is at quasi- 
equilibrium, and in particular, that wind stresses are steady. So the winds used in both L1 and L2 
models, consisting of a seven-day running average of the Navy Corrected Geostrophic Winds 
(Rhodes et al., 1986 and 1989), are sampled every three days with linear interpolation used in time 
between samples. As in the Year 3 simulations, heat fluxes are represented by relaxation to a 
monthly density climatology. 

Experiment 212/83.1 has been selected as the best two-layer simulation (L2). Model 
parameters are: 

° Upper-layer inflow transport = 26 x 106 m3 sec-! (26 Sv), 

° Lower-layer inflow transport = 5 x 106 m3 sec-! ( 5 Sv), 

. Wind stress = seven day Navy Winds, 

° Heat flux = relaxation to monthly density climatology, 

. Horizontal eddy viscosity, A = 200 m2/sec, 

. Horizontal density diffusivity, ARHO = 1000 m2/sec, 

. Grid spacing = 20 by 22 km (0.2 by 0.2 degrees), 

. Upper-layer reference thickness, Hl = 200 m, 

. Lower-layer reference thickness, H2 = 3,450 m, 

. Minimum depth of bottom topography = 300 m, 

° Beta, df/dy = 2*10-11 m-! sec-1, 

. Coriolis parameter at the southern boundary, f = 4.5 x 10-5 sec-!, 

. Gravitational acceleration, g = 9.8 m/sec?2, 

. Reference reduced gravity, g' = .025 m/sec?, 

° Interfacial stress = 0, 

° Coefficient of quadratic bottom stress = .002, and 

° Time step = 0.5 hours. 


The simulation was spun up from rest for 10 years to reach statistical equilibrium and then 
run for another 10 years. All results are from the last 10 years. Figures 18-25 show free surface 
deviation and upper-layer density deviation every 90 days for a typical Loop Current eddy cycle. 
This sequence was forced by winds from 1971 and 1972, but the simulation is not intended to be a 
hindcast, i.e., the circulation is typical of the Gulf, but does not represent the actual Gulf in 1971 
and 1972. The path of the Loop Current eddy, into the southwest Gulf and then north to the Texas 
coast, is typical of large eddies. The simulation includes smaller Loop Current eddies that take a 
path due west, rather than southwest. Figures 26 and 27 show free surface deviation every 90 
days for a typical cycle that includes the reabsorption of an eddy by the Loop Current. Over the 10- 
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Figure 18. Instantaneous view of (a) the free surface deviation and (b) the upper-layer density 
deviation from Experiment 212/83.1 on wind Date 093/1971. 
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Figure 19. Instantaneous view of (a) the free surface deviation and (b) the upper-layer density 


deviation from Experiment 212/83.1 on wind Date 183/1971. 
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Figure 20. Instantaneous view of (a) the free surface deviation and (b) the upper-layer density 
deviation from Experiment 212/83.1 on wind Date 273/1971. 
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Figure 21. Instantaneous view of (a) the free surface deviation and (b) the upper-layer density 
deviation from Experiment 212/83.1 on wind Date 004/1972. 
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Figure 22. Instantaneous view of (a) the free surface deviation and (b) the upper-layer density 
deviation from Experiment 212/83.1 on wind Date 094/1972. 
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Figure 23. Instantaneous view of (a) the free surface deviation and (b) the upper-layer density 
deviation from Experiment 212/83.1 on wind Date 184/1972. 
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Figure 24. Instantaneous view of (a) the free surface deviation and (b) the upper-layer density 
deviation from Experiment 212/83.1 on wind Date 274/1972. 
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Figure 25. Instantaneous view of (a) the free surface deviation and (b) the upper-layer density 
deviation from Experiment 212/83.1 on wind Date 004/1973. 
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Figure 26. Instantaneous view of the free surface deviation from Experiment 212/83.1 on wind 
Dates (a) 005/1974 and (b) 095/1974. 
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Figure 27. Instantaneous view of the free surface deviation from Experiment 212/83.1 on wind 
Dates (a) 185/1974 and (b) 275/1974. 
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year period, six eddies are shed from the Loop Current. The average Loop Current eddy shedding 
cycle in the simulation is perhaps as much as twice that of the actual Gulf. This may be partly due 
to the additional friction from the shelf mixing. The effective distance between the Yucatan and 
Florida Straits is also reduced by the shelf mixing and has also been shown to have the effect of 
increasing the shedding period in the case of highly idealized Gulf of Mexico geometries (Hurlburt 
and Thompson, 1980). 


Figure 28 shows interface deviation and upper-layer density deviation annual climatologies 
over the 10 years for Experiment 121/83.1. These fields are comparable to the depth of the 
thermocline and the temperature at that depth (i.e., Robinson, 1973). Figures 29-32 show 
monthly climatologies for January, April, July, and October. There is a clear seasonal cycle in the 
density climatology, but a 10-year sample length is too short for any significance to be placed in 
most of the differences in the interface deviation fields. 


Experiment 212/80.0 has been selected as the best one-layer simulation (L1). Model 


° Upper-layer inflow transport = 10 x 106 m3 sec-! (10 Sv), 
. Wind stress = seven-day Navy winds, 

. Horizontal eddy viscosity, A = 200 m2/sec, 

. Grid spacing = 20 by 22 km (0.2 by 0.2 degrees), 

. Upper-layer reference thickness, H1 = 3,650 m, 

. Minimum depth of bottom topography = 30 m, 

. Beta, df/dy = 2*10-1! m-! sec-l, 

. Coriolis parameter at the southern boundary, f = 4.5 x 10-5 sec-!, 
. Gravitational acceleration, g = 9.8 m/sec?, 

. Coefficient of quadratic bottom stress = .002, and 

. Time step = 0.5 hours. 


The simulation was spun up from rest for six years to reach statistical equilibrium and then 
run for another 10 years. All results are from the last 10 years. The forcing fields are identical to 
those used in Experiment 121/82 * final 10 years of the two simulations have been merged 
as described above (L2+L1). Fi,..¢ © | shows upper-iayer currents from both experiments (L2 
and L1) on the same day for the Texas shelf region. The 100 and 200 raeter topography depth 
contours are included to show the blending region. Figure 34 shows the corresponding merged 
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Figure 28. Annual cli 


, from 10 years of Experiment 212/83.1 (a) interface deviation and 


limatology 
(b) upper-layer density deviation. 
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Figure 29. Monthly climatology for January from 10 years of Experiment 212/83.1 (a) interface 


deviation and (b) upper-layer density 
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Figure 30. Monthly climatology for April from 10 years of Experiment 212/83.1 (a) interface 
deviation and (b) upper-layer density deviation. 
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Figure 31. Monthly climatology for July from 10 years of Experiment 212/83.1 (a) interface 
deviation and (b) upper-layer density deviation. 
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Figure 32. Monthly climatology for October from: i0 years of Experiment 212/83.1 (a) interface 
deviation and (b) upper-layer density deviation. 
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Figure 33. Upper-layer current vectors and bottom topography contours for the Texas shelf on 
wind Day 001/1969 for (a) one-layer Experiment 212/80.0 and (b) two-layer 
Experiment 212/83. 1. 


CURRENT AND B. TOPOCRAPHY TEXAS SHELF 212321 800 
“AYER = | 12 DATE = S01, 1969 
3iN - : 


: 
‘ j 
30N | , } 
: ———- + © 
' ; —- ase Mceegingihy Maa _ 
; . Oe eet: 
29N “s**e ee ee ee oe 
oo eee werweueeeuteneens 
a o— Pero 
, -?_ ae - 
| 7 , ee 
La Z =< ¢ 
| y “gfe @<enartr*t Cere 4<46444¢€¢¢€¢6€*" © Ene revere TT P4aenet 
fy) $Psacrrrr er eee 2£¢feer*? RP © Pee eee 7 AR ee 
fe ATP Perr ree eer rerere , _ - + * ° > FFP ee aana 
“7+ »>er +e >see > ~€enaaa 
27N ' entt Fe ae ee eee ee © © et tne & patti ecoo ee 
‘ wet? . “7 Ff 7,7 eee wee tee eeaeVecee 
e* a¢4erv77 0? © .* ere ,4aeangeeerrerer ? ° . -_ * #2 
{ 4* *¢>, ? sec A ae al 
i \ s* e*ege *ee - *= « FF > * Ff = > -— © +e. oe 
26N <* *- ** ees ~ 
Ne set hae > = = = — -— © & © © geet e ee eres 
: ge eet 45 ded . 2 ee ee 
4% * & 4apee + > > Oe EHH Ce mee Cees 
} » fe es a paeteeeanes ECE « 
25N ~ FPP 6 4aene 94448 ee ee BRR BORA 
ded ev. 7? 7 4aanaee aah @ 
| ¢ @ @ @ 44eee ** er annanne F © © & BPePPR 44 oe 
0006 a8 vane * tae? Prrveeest VU ves uge FIGURE 33a 
4a ‘ . "<4 
oat ? v,eaaaee th tadadctdrtedadedter *“*f**teers sa? 


MAX SPEED «= 0243 M/S 
NORDA 323 11 -O€C -@8 


CURRENT AND B. TOPOGRAPHY TEXAS SHELF 21242:2 83.1 


LAYER «= | a2” DATE = 001/1969 
JIN ” a T T : ¥ a ' 1 
' T 
] * ‘ I] 
‘ ® od ; 
: | 
30N - - - 
ry 7 a 
i ’ : : ; 
: : --* . oe. 
‘ ‘ _ee Seeeee = : 
; - ‘ 7 eee Ge e@he e@ 
29N “eee wecee Perea tsar eehreccecenae tee j 
° : CN 16 Se om 
‘ °° ©4424 ¢ © ©4706 wae a“ an 


‘ 
Sk akiTEE se, 
+ mere + é 
Faeteee Sy, 


“sereteuee 


td deed | ‘| 
pene teds a a * Se -f 

«¢ 
ae ° 


Hikes 
MH \ 
~e ett ey 


ere eee es 
Cae hte eeeans 
oo ee esccane s 47 
OC at eter eee ant: / r, 


: 3: set Reanaaes Me (i 


go feSSanae ** 
oe eennneer 


yr aee® 


4, 


7] 1 Pppeee 


Z Le eee? 
ih Vif ey 


wifi" 


FIGURE 33b 


mo@oe 323 11 -OfC~-O8 


48 BEST COPY AVAILABLE 


Figure 34. Merged current vectors and bottom topography contours for the Texas shelf on wind 
Day 001/1969: (a) “geostrophic” surface currents and (b) surface currents. 
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field (L2+L1), labeled as "geostrophic” surface current, and the surface current obtained after 
applying the perturbation analysis to obtain a detailed vertical profile (L2+L1+P). Figures 35-37 
show the same merged fields at 90-day time intervals to illustrate the seasonal cycle on the shelf. 
Figures 38 and 39 are one year later than Figures 34 and 35 respectively, illustrating interannual 
variability. 

The present Oil Spill Risk Analysis system at Minerals Management Service calculates the 
surface current due to local wind forcing by taking 3.5% of the wind speed at an angle to the right 
of the wind of between 0 and 30 degrees depending on wind speed (Samuels et al., 1982). This 
current is then added to the geostrophic component to obtain the total surface current. Figure 40 
shows the surface current for a day in January calculated via the perturbation analysis and by the 
3.5% rule (L2+L1+P vs. L2+L1+3.5%). Figures 41 and 42 show the corresponding wind 
stresses and the surface currents minus the geostrophic component (P only versus 3.5% only). 
Figures 43-45 show the same fields for a day in April. The perturbation analysis typically 
generates a current component that is about 45 degrees to the right of the local wind, which is 
further to the right than the 3.5% rule. This is consistent with the Ekman spiral produced by the 
analytic analysis and is appropriate for a steady wind. The different angles, and the fact that the 
perturbation analysis gives smaller magnitudes, is partly due to the absence of the Stoke's drift 
current component from the analysis. The perturbation analysis is also sensitive to the drag 
coefficient used to convert wind speeds to stresses. In this case, the drag coefficient is 0.00012 
[see Rhodes et al. (1989)]. These winds were not generated as part of this contract; if they were, it 
is likely that greater attention would have been paid to what drag coefficient to use. The 
geostrophic surface currents (L2+L1) and the wind stresses have been made available, in addition 
to the surface currents from the perturbation analysis (L2+L1+P). So the 3.5% rule, or some other 
alternative technique, can be used to obtain the currents due to local wind effects. Figures 46-52 
cover the same period as Figures 33-39, but for the Florida shelf. Obviously, this brief 
comparison is not sufficient to decide on the relative merits of the two approaches. In any case, 
both are very gross approximations to the actual effect of winds on the surface current. 


Figure 53 shows the annual surface current climatology for the full 10-year simulation 
(L2+L1+P) and a corresponding climatology based on historical ship drift data (N) provided by the 
Naval Oceanographic Office. The observational climatology bins data into one-degree squares, so 
the full resolution of this data set is shown. The model currents have been subsampled in space to 
improve plot clarity. Given the wide variation in horizontal resolution between these two data sets, 


Figure 35. Merged topograp 
Day 091/1969: (a) “geostrophic” surface currents and (b) surface currents. 
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Figure 36. 


Merged current vectors and bottom 


hy contours for the Texas shelf on wind 


: topograp 
Day 181/1969: (a) "geostrophic” surface currents and (b) surface currents. 
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Figure 37. Merged current vectors and bottom 


y contours for the Texas shelf on wind 


Day 271/1969: (a) “geostrophic” surface currents and (b) surface currents. 
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Figure 38. Merged current vectors and bottom topography contours for the Texas shelf on wind 
Day 002/1970: (a) “geostrophic” surface currents and (b) surface currents. 
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Figure 39. Merged current vectors and bottom topography contours for the Texas shelf on wind 
Day 092/1970: (a) “geostrophic” surface currents and (b) surface currents. 
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Figure 40. Merged current vectors and bottom topography contours for the Texas shelf on wind 
Day 001/1969: (a) surface currents and (b) surface currents using the 3.5% rule. 
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Figure 41. a a 
1/1969. 
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Figure 42. (a) Wind stresses and (b) local wind-induced currents from the 3.5% rule for the 


Texas shelf on wind Day 001/1969. 
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Figure 43. Merged current vectors and bottom topography contours for the Texas shelf on wind 
Day 091/1969: (a) surface currents and (b) surface currents using the 3.5% rule. 
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Figure 44. (a) Wind stresses and (b) local wind-induced currents for the Texas shelf on wind Day 
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Figure 45. (a) Wind stresses and (b) local wind-induced currents from the 3.5% rule for the 
Texas shelf on wind Day 091/1969. 
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Figure 46. Upper layer current vectors and bottom topography contours for the Florida shelf on 
wind Day 001/1969 for: (a) one-layer Experiment 212/80.0 and (b) two-layer 
Experiment 212/83.1. 
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Figure 47. Merged current vectors and bottom topography contours for the Florida shelf on wind 
Day 001/1969: (a) “geostrophic” surface currents and (b) surface currents. 
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Figure 48. Merged current vectors and bottom topography contours for the Florida shelf on wind 
Day 091/1969: (a) “geostrophic” surface currents and (b) surface currents. 
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Figure 49. Merged current vectors and bottom topography contours for the Florida shelf on wind 
Day 181/1969: (a) "geostrophic” surface currents and (b) surface currents. 
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Figure 50. Merged current vectors and bottom topography contours for the Florida shelf on wind 
Day 271/1969: (a) “geostrophic” surface currents and (b) surface currents. 
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Figure 51. Merged current vectors and bottom topography contours for the Florida shelf on wind 
Day 002/1970: (a) “geostrophic” surface currents and (b) surface currents. 
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Figure 52. Merged current vectors and bottom topography contours for the Florida shelf on wind 
Day 092/1970: (a) “geostrophic” surface currents and (b) surface currents. 
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Figure 53. Annual surface current climatology (a) from the 10-year merged simulation and (b) 
from historical observations. 
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a side-by-side comparison is more appropriate than a plot of the two fields differenced together. 
Both fields include a mean Loop Current, but it bends further westward in the simulation. Both 
show westward zonal flow at 22N to 23N bending northward near the Mexico/Texas coast. These 
Currents are stronger in the observations, but the simulation only includes seven-day wind stress 
averages and the simulated mean flow might be stronger with 12 hourly winds. Figure 54 
compares the annual surface current and annual geostrophic surface current climatologies. The 
zonal flow at 22N is weaker in the geostrophic case, suggesting this is largely a wind driven 
response. Figures 55-58 show seasonal climatologies. In all seasons the simulation shows a 
cyclonic gyre in the far soutiiwest Gulf. This feature is not seen in the observations; however, it 
would be barely resolvable on a one-degree grid and there are probably relatively few observations 
in the region. On the other hand, the wind stresses used by the model are also open to question in 
that region (Rhodes et al., 1989). Figures 59-62 show a representative instantaneous field and its 
difference from the seasonal mean for each season. The differences are taken with respect to the 
model mean, rather than the mean from observations because of the wide variation in horizontal 
resolution between these two data sets. The size of the anomaly is comparable to that of the 
original field particularly in deep water. 

The comparisons, begun in Year 2 (Wallcraft, 1986), between statistical measures of the 
simulated currents and those from actual field measurements have been continued. The ocean 
model simulation is not a hindcast (i.¢., it does not claim to represent the actual situation in the 
Gulf over any given time period), so direct day-by-day comparison with observations would be 
inappropriate. This is unfortunate since a direct comparison would provide far more information 
on the quality of the simulation than the statistical measures used here. Figure 63 shows the 
locations of moored current meter arrays deployed as part of MMS's Gulf of Mexico Physical 
Oceanography Program (SAIC, 1986). Very good data has been obtained from these moorings, 
but statistical comparisons are made more difficult by the fact thas the data record at many moorings 
extends for no longer than one year. The average Loop Curren’ eddy shedding period is about one 
year, but it can vary from 6-18 months. So a one-year data record is too short to obtain reliable 
statistics about Loop Current related phenomena, which means in practice that the only place where 
the data record is likely to be long enough for good long-term statistics is over the shelf. Figure 64 
shows the full 10-year velocity component spectra from the simulation at mooring MOF, which is 
on the inner shelf so the one-layer model (L1) results apply. There is very little energy below 10 
days, as would be expected from seven-day averaged wind stress forcing, and there is a strong 
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Figure 54. Annual surface current climatology (a) from the 10-year merged simulation and (b) 
" climatology from the merged simulation. 
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Figure 55. Winter surface current climatology (a) from the 10-year merged simulation and (b) 
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Figure 56. 


Spring surface current climatology (a) from the 10-year merged simulation and (b) 


from historical observations. 
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Figure 57. Summer surface current climatology (a) from the 10-year merged simulation and (b) 
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Figure 58. Fall surface current cli a) from the 10- 
historical it climatology (a) year merged simulation and (b) from 
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Figure 59. On wind Day 001/1969 (a) surface current from the merged simulation and (b) surface 
Current anomaly versus the winter climatology. 
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Figure 60. On wind Day 091/1969 (a) surface current from the merged simulation and (b) surface 
current anomaly versus the spring climatology. 
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Figure 61. On wind Day 181/1969 (a) surface current from the merged simulation and (b) surface 
Current anomaly versus the summer. 
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Figure 62. On wind Day 271/1969 (a) surface current from the merged simulation and (b) surface 
Current anomaly versus the fall climatology. 
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Figure 63. Map of the Gulf of Mexico showing the locations of current meter arrays from the 
Gulf of Mexico Physical Oceanography Program. 
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Figure 64. Velocity component spectra for 10 years from the simulation at location MOF (a) u 
component and (b) v component. 
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annual signal in the cross shelf, but not the the long shelf direction. In contrast to the case studied 
in the second year (Wallcraft, 1986) comparison (Figure 11), there is no peak at about 100 days in 
the latest simulation. This is additional evidence that the original signal was due to Loop Current 
effects that cannot be present in the merged simulation at this location. Figures 65 and 66 show 
velocity component spectra at location MOD for the observations and the simulation. As in all 
cases, there is very little energy below seven days in the simulation. Both show a peak at about 
50-100 days, but the simulation shows a strong peak at 20 days that is not present in the 
observations. Figures 67 and 68 are for location MOE that has one of the longest set of 
observations (1,090 days). The observations and simulation are in good agreement although the 
10-year simulation shows energy at multiyear periods that is not present in the observations. 

Figures 69 and 70 are for the simulation only at locations MOA and MOC. They show a 
strong peak at about two years, which is the signature of the Loop Current cycle in this simulation 
(six eddies shed in 10 years). Figures 71 and 72 are again for the simulation only, but this time for 
locations MFF and MGG in the central Gulf. A strong Loop Current eddy signature can be seen 
here too. 


Wind stress, surface current (L2+L1+P), and geostrophic surface current (L2+L1) fields 
from the merged simulation (sampled every three days for 10,750 and 1,600 meter depths, 
sampled every six days for the same 10-year period) have also been delivered, as have detailed 
vertical current profiles (L2+L.1+P) at the locations of all the moored buoys used in Years 1-5 of 
MMS's recently concluded Gulf of Mexico Physical Oceanography Program. Figures 73-193 
show all the surface current fields (L2+L1+P) delivered for one full Loop Current eddy cycle, 
which is about 20% of the surface current fields in the full 10 years. 


Vv. CONCLUSIONS 


This project was designed from the outset to be modest in scope. The limitations of the 
original layered model were known from the beginning (JAYCOR, 1983), but the fact that the 
shelf circulation is vital to the success of any simulation used in oil spill risk analysis was not fully 
appreciated at first. Attempts to improve the multilayered model's skill over the shelf, while 
maintaining the project's modes f:nding level, have proved largely unsuccessful. The final 
product does not have several of the major enhancements to the model that were originally 
proposed; and instead, uses two related models with the same forcing, but with no other coupling 
between them for the shelf and deep waiter circulation. If a more ambitious project were to choose 


Figure 65. Velocity component spectra for 362 days from the obs=rvations at location MOD (a) u 
component at 60 m and (b) v component at 60 m. 
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Figure 66. Velocity component spectra for one year from the simulation at location MOD (a) u 


component and (b) v component. 
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Figure 67. Velocity component spectra for 1090 days from the observations at location MOE (a) u 
component at 100 m and (b) v component at 100 m. 
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Figure 68. Velocity component spectra for 10 years from the simulation at location MOE (a) u 
component and (b) v component. 
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Figure 69. Velocity component spectra for 10 years from the simulation at location MOA (a) u 


component and (b) v component. 
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Figure 70. Velocity component spectra for 10 years from the simulation at location MOC (a) u 
component and (b) v component. 
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Figure 71. Velocity component spectra for 10 years from the simulation at location MFF (a) u 
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Figure 72. Velocity component spectra for 10 years from the simulation at location MGG (a) u 


component and (b) v component. 
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Figures 73-193. Surface currents from the merged simulation every three days for a full Loop 
Current eddy cycle. The wind day for each figure is given by the "DATE" in the 
top right-hand comer of each plot. 
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to use a layered model in the Gulf, it should probably look first to flux-corrected transport 
techniques for producing improved simulations; however, such a project would also have the 
option of using several other kinds of ocean models that did not use layers in the vertical. 


The final geostrophic surface currents (L2+L1) and surface currents (L2+L1+P) from this 
Study provide a vastly improved realization of Gulf circulation than the existing alternatives; 
however, there are several limitations in the system design that suggest better simulations may be 
possible, such as the following approaches: 


° The most obvious limitation of the deep water circulation is that the simulation's 
typical Loop Current eddy shedding period, at about two years, is twice as long as 
the generally-accepted average period for the Gulf. Other two-layer simulations 
have had shorter eddy shedding periods, but were less suitable for merging with the 
one-layer results. 


° The transport through the Yucatan Straits was held constant in both magnitude and 
profile throughout all the simulations. The only possible source for a detailed long- 
term time series of currents through the Straits would be a numerical model of the 
entire North Atlantic forced by observed winds. The effects of variations in flow 
through the Yucatan Straits on Gulf circulation are not known; if they are indeed 
2 ey eo 
the ; 


° A better simulation of shelf circulation could almost certainly be obtained by using a 
model designed specifically for that purpose. 
° The validity of blending two simulations at all is open to question. 


The most important model-related question for future field measurement programs is the 
nature and the degree of interaction between on- and off-shelf circulation. It might be that such a 
program would demonstrate that on- and off-shelf circulation can indeed be treated separately 
without significantly degrading simulation skill. At the very least it would quantify the kinds of 
interactions a model, or suite of models, must have to produce better simulations. 


It is probably not possible to produce long-term simulations of the deep water circulation in 
any section of Gulf of Mexico without modeling the entire Gulf at a horizontal resolution of at least 
25 km. Too little is known about the interactions between on- and off-shelf circulation to be 
confident that simulations that decouple these regimes can, even in principle, capture all the 
important variations in the current field. Therefore, the ideal long-term simulation in the Gulf 
would come from a model of the entire Gulf of Mexico or perhaps a system of tightly coupled 
models of the Gulf that demonstrated simulation skill in both the deep water and on the shelf. The 
problem with this approach is that the vertical scales and, to a lesser extent, the horizontal scales 
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required to model on- and off-shelf processes are entirely different. It is extremely difficult to 
include both sets of scales and the transition between them in one model at the same time. There 
are no models that have demonstrated the ability to do this successfully in any comparable oceanic 
region and there are only a few model designs that might, theoretically, be able to so given 
sufficient resources. 


On the horizon, there is an alternative to simulations of any kind. In the atmosphere the 
best source of wind fields are the prediction centers that produce atmospheric forecasts every day. 
Any long-term forecast sequences will provide better statistics than a long-term numerical 
simulation of the atmosphere. Similarly a long-term sequence of ocean forecasts would be 
superior to any simulation of the Gulf. The routine forecasting of ocean currents on the scales 
necessary for good statistics in the Gulf of Mexico should be a reality by 1994. These forecasts 
will immediately provide a valuable source of verification data for numerical simulations and will 
entirely supplant the need for simulations within 10 years of their introduction. 
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